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0.2 Zusammenfassung 



Im Minimal Supersymmetrischen Standardmodell (MSSM) mit komplexen Parametern 
ergeben Einschleifen-Strahlungskorrekturen des Zerfalls eines Stop in ein Bottom-Quark 
und ein Chargino eine CP verletzende Asymmetric der Zerfallsbreite 

Wir fiihren einc detaillierte, numerische Analyse von 6^^ als auch von 5*"^ x BR (wobei 
BR das Verzweigungsvcrhaltnis dcs Zerfalls ist) durcli, wobei wir die Abhangigkeit von 
den beteiligten Parametern und komplexen Phasen untersuchen. Dabei nehmen wir die 
Yukawa-Kopplungskonstanten des Top- und Bottom-Quarks als laufend an. Wir beriick- 
sichtigen die Einschrankungen der Parameter, welche von der experimentcllcn Obergren- 
ze des elektrischen Dipolmoments des Elektrons ausgehen, indem wir das Dipolmoment 
wahrcnd dor Analyse automatisch mitbcrcchncn und damit kontrollicrcn. 
Als Rcsultatc crlialten wir fiir die Asymmetric 6'^"^' einen Wert bis zu ~ 24%, abliangig 
vom gewahlten Punkt im Parameterraum. Die kombinierte Grofie x BR wird bis zu 
~ 3.5 % groB. 

Wir kommentieren auch die Moglichkcit ciner Messung dicser Asymmetric am Large 
Hadron Collider (LHC) dcs CERN. Es wird moglich sein, unsere Asymmetric der Zer- 
fallsbreite 6'"^ am LHC zu messen. 
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0.3 Abstract 



In the Minimal Supersymmetric Standard Model (MSSM) with complex parameters, 
one-loop corrections to the decay of a stop into a bottom-quark and a chargino lead to 
the CP violating decay rate asymmetry 

T^(i,^bxt)-T-{i*^bxr) 

We perform a detailed numerical analysis of 6'-''^ and also 6'-''^ x BR (where BR is the 

branching ratio of the decay) for ti ^ bxt and t2 — > bxt, analyzing the dependence 

on the parameters and complex phases involved. In addition, we take the Yukawa 

couplings of the top- and bottom-quark running. We account for the constraints on 

the parameters coming from the experimental limit of the electric dipole moment of the 

electron by calculating and thus checking it automatically along the way. 

We obtain as results that the asymmetry d^^ rises up to ~ 24 %, depending on the point 

in parameter space. The combined quantity 6^^ x BR reaches up to ~ 3.5 %. 

We also comment on the feasibility of measuring this asymmetry at the Large Hadron 

Collider (LHC) at CERN. It will be possible to measure our decay rate asymmetry S'^^ 

at LHC. 
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Chapter 1 
Introduction 



Although the Standard Model (SM) of elementary particle physics is a remarkably suc- 
cessful description of presently known phenomena, it has to be extended to describe 
physics at high energies properly. Among the theoretical and phenomenological issues 
that the SM fails to address properly is the baryon asymmetry of the universe [Ij. 
Phenomenologically, there are many reasons to believe that we live in a baryon asym- 
metric universe, so that there exists much more matter than anti-matter. One strong 
piece of evidence comes from the acoustic peaks — early universe baryon-photon plasma 
oscillations — inferred from Cosmic Microwave Background measurements (see e.g. ^), 
which give the baryon-to-photon ratio 



in which s is the entropy density (roughly the photon density), and Ub and are the 
number densities of baryons and antibaryons, respectively. This data agrees well with 
big bang nucleosynthesis, which requires the baryon-to entropy density ratio to be (see 



The problem of baryogenesis [5] is to explain the origin of this small number starting 
from the natural initial condition of 77 = 0, which in most cases is attained at high 
enough temperatures. 

Assuming CPT is preserved, there are three necessary conditions for baryogenesis, usu- 
ally referred to as the Sakharov requirements [5]: 

1. Baryon number violation 

2. Departure from thermal equilibrium 

3. Charge (C) and Charge-Parity (CP) violation. 

Although the Sakharov criteria can be met in the SM, the baryon asymmetry generated 
at the electroweak phase transition is too small. However, supersymmetric extensions 
of the SM can contain new sources of CP violation which can increase and thus might 
explain the baryon asymmetry of the universe. 
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CHAPTER 1. INTRODUCTION 



Supersymmetric models introduce superpartners to every particle, that just differ in 
the spin quantum number. Since we do not observe such superpartners to the known 
particles with the same mass, supersymmetry must be broken. This supersymmetry 
breaking leads to a couple of new parameters. If these parameters are chosen to be 
complex, supersymmetric models contain new sources of CP violation. In the Minimal 
Supersymmetric Standard Model (MSSM) — the most promising extension of the SM 
— the U (1) and SU (2) gaugino (superpartner of gauge bosons) mass parameters Mi and 
M2, respectively, the higgsino (superpartner of Higgs-boson) mass parameter /x, as well 
as the trilinear couplings Af (corresponding to a fermion f) may be complex. Radiative 
corrections at one-loop level can then lead to new CP violating asymmetries in addition 
to the small CP asymmetries in the SM coming from the Cabibbo-Kobayashi-Maskawa- 
(CKM)-matrix. 

Although large complex phases are desirable to explain baryogenesis [6], [TJ [8|, [9] , there are 
constraints on these phases coming from the experimental limits on the electric dipole 
moments (EDMs) of the electron, neutron and Hg. Especially the complex phase of fi 
is highly constrained for a typical SUSY mass scale of the order of a few hundred GeV. 
In this thesis we study the CP violating decay rate asymmetry 

gcp ^ T^{i.^bxt)'T'{i:^bxr) 3. 

of the decay of a stop (bosonic superpartner of the top-quark) into a bottom-quark and 
a chargino (mass eigenstate of the superpartners of charged Higgs bosons and charged 
gauge bosons) in the MSSM with complex parameters at full one-loop level. 
The asymmetry is of course zero if CP is conserved and also vanishes at tree-level in the 
case of CP violation. The complete list of all graphs at one-loop level which contribute to 
this asymmetry can be found in the Appendix \M They give a contribution only if they 
have an absorptive part, i.e. at least a second decay channel of ti must be kinematically 
possible in addition to that into • spite of the overall 47 graphs one, however, 
expects that two graphs will dominate. They are the two possible processes involving 
a gluino g (see Chapter [H]). Because the gluino couples like its superpartner the gluon 
with the strong interaction force, these contributions are expected to dominate over all 
others, if the decay channel ti ^ tg is open (m^-. > nit + ^g)- 

As a loop-level quantity the decay rate asymmetry 5'"'^ depends on the phases of all 
complex parameters involved. One, however, expects that the dependence on the phase 
of At and A^, is strongest (taking /x real because of the stringent EDM constraints). We 
neglect the mixing in the neutral Higgs sector due to the complex trilinear couplings of 
the third generation {At, Ab, Ar) (see e.g. [TOl [HI [12]) because they are very small [13]. 

We perform a detailed numerical analysis for ti bxt and t2 ^ bxt analyzing the 
dependence on the parameters and phases involved. In addition, we take the Yukawa 
couplings of the top- and bottom-quark running. We account for the constraints coming 
from the EDM of the electron by calculating and thus checking it automatically along 
the way. We also comment on the feasibility of measuring this asymmetry at the Large 
Hadron Collider (LHC) at CERN. 
The work is organized as follows: 
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• In Chapter 2 we give an introduction into the basic concept of supersymmetry 
(SUSY). 

• In Chapter 3 we derive supersymmetric Lagrangians as well as the soft supersym- 
metry breaking interactions. 

• In Chapter 4 we present the particle content of the Minimal Supersymmetric Stan- 
dard Model (MSSM). 

• In Chapter 5 the mass matrices of the sfermions, charginos and neutralinos are 
derived. 

• In Chapter 6 we list the relevant couplings of the most important interactions 
involved. 

• Chapter 7 contains some definitions which are important for the study of CP 
violating asymmetries. 

• Chapter 8 provides a deeper insight into CP violation and the decay rate asym- 
metry S'"^. 

• In Chapter 9 we list the most important processes who are expected to yield the 
highest S'^^. 

• In Chapter 10 we discuss the numerical results and present the conclusions. 

• Appendix A holds the complete list of all processes at full one-loop level who can 
contribute to 5^^. 

• In Appendix B we define the Passarino-Veltman integrals and provide them with 
a special argument set. 

• In Appendix C and D we derive a generic structure and a tree-level coupling, 
respectively. 

• Appendix E shows how to transform from Weyl to Dirac spinors and vice versa. 

• Finally, in Appendix F we derive the electric dipole moment (EDM) of a fermion 
and calculate the EDM of the electron. 



CHAPTER 1. INTRODUCTION 



Chapter 2 

Supersymmetry (SUSY) 



2.1 Beyond the Standard Model 

The Standard Model of elementary particle physics (SM) [m |15l [16] is a spectacularly 
successful theory of the known particles and their electroweak and strong forces. The SM 
is a gauge theory, in which the gauge group SU{3)c x SU{2)l x U{1)y is spontaneously 
broken to SU{3)c x U{1)em by the non- vanishing vacuum expectation value (VEV) of 
a fundamental scalar field, the Higgs field, at energies of order 100 GeV. Although the 
SM provides a correct description of all known microphysical non-gravitational phenom- 
ena (except neutrino masses and oscillations), there are a number of theoretical and 
phenomenological issues that the SM fails to address adequately [1]: 

• Hierarchy problem. Phenomenologically the mass of the Higgs boson associated 
with electroweak symmetry breaking must be in the range of (9(100 GeV). How- 
ever, radiative corrections to the Higgs mass are quadratically dependent on the 
ultra-violet cutoff A, since the masses of fundamental scalar fields are not protected 
by chiral or gauge symmetries. The "natural" value of the Higgs mass is therefore 
of 0{A) rather than (9(100 GeV), leading to a destabilization of the Higgs mass 
and the hierarchy of the mass scales in the SM. To achieve m ~ (9(100 GeV) it is 
necessary to fine-tune the scalar mass-squared parameter rriQ ~ A^ of the funda- 
mental ultraviolet theory to a precision of m^/A^. If, for example, A = 10^^ GeV 
and m = 100 GeV, the precision of fine-tuning must be 10~^^, which is very un- 
natural. 

• Electroweak symmetry breaking. In the SM, electroweak symmetry breaking is 
parameterized by the Higgs boson h and its general potential F = /i^|/ip + A|/;.|^. 
In order to make this symmetry breaking happen, one has to set /i^ < by hand. 
This postulate is rather artificial. 

• Gauge coupling unification. The idea that the gauge couplings undergo renormal- 
ization group evolution in such a way that they meet at a point at a high scale 
lends credence to the picture of grand unified theories (GUTs) and certain string 
theories. However, precise measurements of the low energy values of the gauge 
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CHAPTER 2. SUPERSYMMETRY (SUSY) 



couplings demonstrated that the SM cannot describe gauge couphng unification 
(see e.g. [17]) accurately enough to imply it is more than an accident. 

• Family structure and fermion masses. The SM does not explain the existence of 
three families and can only parameterize the strongly hierarchical values of the 
fermion masses. Massive neutrinos imply that the theory has to be extended, 
as in the SM the neutrinos are strictly left-handed and massless. Right-handed 
neutrinos can be added, but achieving ultra-light neutrino masses from the seesaw 
mechanism [T8l [T9] requires the introduction of a new scale much larger than 
0(100 GeV). 

• Cosmological challenges. Several difficulties are encountered when trying to build 
cosmological models based solely on the SM particle content. As already mentioned 
in the introduction, the SM cannot explain the baryon asymmetry of the universe. 
The SM also does not have a viable candidate for the cold dark matter of the 
universe, nor a viable infiaton. The most difficult problem the SM has when 
trying to connect with the gravitational sector is the absence of the expected scale 
of the cosmological constant. 

Therefore, the Standard Model must be extended to be valid at higher energies. Theories 
with low energy supersymmetry have emerged as the strongest candidates for physics be- 
yond the SM. The main idea behind supersymmetry (SUSY) is an underlying sjTnmetry 
between bosons and fermions. In the simplest supersymmetric world, each SM particle 
has a superpartner which differs only in the spin by 1/2 and is related to the original 
particle by a supersymmetry transformation. 

Since no superpartners were discovered yet, SUSY is not an exact symmetry. They seem 
to have a significant higher mass, which can be explained by spontaneous breaking of 
the supersymmetry. In other words, the underlying model should have a Lagrangian 
density that is invariant under supersymmetry, but a vacuum state that is not. 

2.2 Predictions and Successes of SUSY 

The main reasons that low energy supersymmetry is taken very seriously are not its 
elegance or its likely theoretical motivations, but its successful explanations and predic- 
tions. Of course, these successes may just be remarkable coincidences because there is 
as yet no direct experimental evidence for SUSY. Superpartners and a light Higgs boson 
must be discovered or demonstrated not to exist at the Large Hadron Collider (LHC) 
at CERN. The main successes are as follows [U |20] : 

• Hierarchy problem. SUSY provides a solution to the hierarchy problem by protect- 
ing the Higgs mass from large radiative corrections coming from heavy particles. 
Due to the symmetry between bosons and fermions, each contribution to the radia- 
tive correction can be completely canceled in the case of exact SUSY. In the case 
of broken SUSY the cancelation is not complete but still strong enough to render 
the corrections to the Higgs mass harmless, if the masses of the new superpartners 
are < C(l TeV). 
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• Radiative electroweak symmetry breaking. With plausible boundary conditions at 
a high scale, low energy supersymmetry can provide the explanation of the origin 
of electroweak symmetry breaking ^211 [221 ISSl [21]. To oversimplify a little, the 
SM effective Higgs potential has the form V = m?h'^ + A/i^. First, supersymmetry 
requires that the quartic coupling A is a function of the f/(l)y and SU(2)c gauge 
couplings A = {g'"^ + 5'^)/2. Second, the m? parameter runs to negative values at 
the electroweak scale, driven by the large top quark Yukawa coupling. Thus the 
"Mexican hat" potential with a minimum away from /i = is derived rather than 
assumed. 

• Gauge coupling unification. In contrast to the SM, the MSSM allows for the 
unification of the gauge couplings, as first pointed out in the context of GUT 
models by [25l [261 EZ]- The extrapolation of the low energy values of the gauge 
couplings using renormalization group equations and the MSSM particle content 
shows that the gauge couplings unify at the scale Mq ^ 3 x 10^^ GeV [281 ISS [301 
[31]. 

• Cold dark matter. In supersymmetric theories, the lightest superpartner (LSP) can 
be stable. This stable superpartner provides a nice cold dark matter candidate [521 
[33]. Simple estimates of its relic density are of the right order of magnitude to 
provide the observed amount. LSPs were noticed as good candidates before the 
need for nonbaryonic cold dark matter was established. 

• Supergravity. Gauged SUSY includes a coupling between gravity and matter. The 
invariance of the Lagrangian density under a local supersymmetry transformation 
leads to a quantized form of Einsteins general relativity. However, like all known 
theories that include general relativity, supergravity is still non-renormalizable as 
a quantum field theory. 

• Baryon asymmetry. As already mentioned in the introduction, SUSY may help 
explain the baryon asymmetry of the universe. At least three different approaches 
can provide the observed baryon asymmetry: (i) generating the asymmetry at 
the electroweak phase transition via the electroweak baryogenesis mechanism (see 
e.g. [HI [8]), (ii) generating it via leptogenesis, and (iii) the Affleck-Dine mecha- 
nism [3l] using the decay of a scalar field in the early universe into matter. All 
three mechanism need the existence of CP violation, as the Sakharov criteria de- 
mand. The CP violation in the SM is too small to explain baryogenesis, but SUSY 
can provide sufficient additional sources of CP violation. 

Supersymmetry has also made several correct predictions [T]: 

• Supersymmetry predicted in the early 1980s that the top quark would be heavy [351 
[55] . because this was a necessary condition for the validity of the electroweak 
symmetry breaking explanation. 

• Supersymmetric grand unified theories (GUTs) with a high fundamental scale 
accurately predicted the present experimental value of sin^ 6y/ before it was mea- 
sured [261 [ig [371 [38]. 
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• Supersymmetry requires a light Higgs boson to exist [391 HO] , consistent with cur- 
rent precision measurements, which suggest Mh < 200 GeV. 

• When the Large Electron-Positron Collider (LEP) began to run in 1989 it was 
recognized that either LEP would discover superpartners if they were very light 
or there would be no significant deviations from the SM (all supersymmetry ef- 
fects at LEP are loop effects and supersymmetry effects decouple as superpartners 
get heavier). In nonsupersymmetric approaches with strong interactions near the 
electroweak scale it was natural to expect significant deviations from the SM at 



Together these successes provide a powerful indirect sign that low energy supersymme- 
try is indeed part of the correct description of nature. 

Remarkably, supersymmetry was not invented to explain any of the above physics. Su- 
persymmetry was discovered as a beautiful theory and was studied for its own sake in 
the early 1970s. Only after several years of studying the theory did it become clear that 
supersymmetry solved the above problems, one by one. Furthermore, all of the above 
successes can be achieved simultaneously, with one consistent form of the theory and its 
parameters. Low energy supersymmetry also has no known incorrect predictions; it is 
not easy to construct a theory that explains and predicts certain phenomena and has 
no conflict with other experimental observations. 



A supersymmetry transformation turns a bosonic state into a fermionic state, and vice 
versa [31]. The operator Q that generates such transformations must be an anticom- 
muting spinor with 



Using the Minimal Supersymmetric Standard Model (MSSM), we only deal with a = 1 
supersymmetry, that means only one generator-pair is needed. Spinors are intrinsically 
complex objects, so (the hermitian conjugate of Q) is also a symmetry generator. 
Because Q and are fermionic operators, they carry spin angular momentum 1/2, so it 
is clear that supersymmetry must be a spacetime symmetry. The possible forms for such 
symmetries in an interacting quantum field theory are highly restricted by the Haag- 
Lopuszanski-Sohnius extension of the Coleman- Mandula theorem |42l|43]. For realistic 
theories like the SM, this theorem implies that the generators Q and must satisfy an 
algebra of anticommutation and commutation relations with the form 



LEP. 



2.3 The Superalgebra 



Q I Boson) = |Fermion) 



(5|Fermion) = | Boson) . 



(2.1) 



{Qa,Q(^} = {QlQl} = 0, 

[Q^,P^] = [QlP>^] = 




(2.3) 
(2.4) 



where P'^ is the four-momentum generator of spacetime translations, the indices a,a = 
1, 2 are the left- and right-handed spinor indices and a stands for the Pauli matrices. 
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The single-particle states of a supersymmetric theory fall into irreducible representations 
of the supersymmetry algebra, called supermultiplets. Each supermultiplet contains 
both fermion and boson states, which are commonly known as superpartners of each 
other. The squared-mass operator commutes with the operators Q, and with all 
spacetime rotation and translation operators, so it follows immediately that particles 
inhabiting the same irreducible supermultiplet must have equal eigenvalues of P^, and 
therefore equal masses. The supersymmetry generators Q, also commute with the 
generators of gauge transformations. Therefore particles in the same supermultiplet 
must also be in the same representation of the gauge group, and so must have the same 
electric charges, weak isospin, and color degrees of freedom. 

Each supermultiplet contains an equal number of fermion and boson degrees of freedom. 
There are in total two different types of supermultiplets possible. The so called chiral 
supermultiplet includes a two-component spin-1/2 Weyl fermion and a complex spin- 
field, whereas the gauge supermultiplet combines a spin-1/2 gaugino (the fermionic 
superpartner of a gauge boson) and a spin-1 gauge boson. 

In a supersymmetric extension of the SM, each of the known fundamental particles is 
thus in either a chiral or a gauge supermultiplet, and must have a superpartner with 
spin differing by 1/2 unit. What remains is to decide exactly how the known particles 
fit into the supermultiplets, and to give them appropriate names. This will be done in 
Chapter HI but first we want to introduce the supersymmetric Lagrangian density and 
the (soft) SUSY breaking terms. 
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Chapter 3 

Supersymmetric Lagrangians 



In this chapter we will describe the construction of supersymmetric Lagrangians and 
soft supersymmetry breaking terms [HI I20j^ . We can then apply these results to the 
special case of the MSSM. 



3.1 A Free Chiral Supermultiplet 

The simplest supersymmetric model we can build is the massless, non interacting Wess- 
Zumino model [H], which describes a single chiral supermultiplet, consisting of a left- 
handed two-component Weyl fermion ip and its superpartner, a complex scalar field 0. 
The simplest action we can write down is 

S d X (/^scalar ~l~ -^fermion) ; ('^■-^) 

Scalar = , (3.2) 

Aermion = llp^a'^d^lp. (3.3) 

A supersymmetry transformation should turn the scalar boson field into something 
involving the fermion field ipa- The simplest possibility for the transformation of the 
scalar field is 

6(l) = eip, 5(p* = , (3.4) 

where is an infinitesimal, anticommuting, two-component Weyl fermion object param- 
eterizing the supersymmetry transformation. As we only discuss global supersymmetry, 
e° is a constant, satisfying d^e°' = 0. The scalar part of the lagrangian now transforms 
as 

^Acalar = C^'^^ ^^0* + C^^^V^t . (3.5) 

We would like for this to be canceled by 5£fermion, at least up to a total derivative, so 
that the action will be invariant under the supersymmetry transformation. This leads 
to only one possibility (up to a multiplicative constant): 

= -t{a^e^)o. d,ct> , = i{ea^)^ d^ct>* . (3.6) 



^In this work we use the signature (+, — , — , — ) of the spacetime metric, contrary to the convention 
found in [411 [20j. On this account some of the terms in our equations have different algebraic signs. 
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With this guess, one immediately obtains 

^-Cfermion = -ea'^d^S^V^ 9^0* + ^^CT^'e^ 9^9,0 . (3.7) 
After applying the Pauli matrix identities 

[^M^- + a'^a^]^'^ = 27]^''6^ , [a^^a" + aV^]^^ = 2r7'^^5f , (3.8) 
and using the fact that partial derivatives commute {d^d,y = 0^0^^), it takes the form 

{ea^a^^/j d^(f)* - 9^0* - eV^ d^cj)) . (3.9) 

The first two terms here just cancel against SCscaiar, while the remaining contribution is 
a total derivative. So we arrive at 

= J d^X (5£scalar + ^Aermion) = , (3.10) 

justifying our guess of the numerical multiplicative factor made in Eq. ( 13. 6p . 
We must also show that the supersymmetry algebra closes; in other words, that the com- 
mutator of two supersymmetry transformations parameterized by two different spinors 
ei and €2 is another symmetry of the theory. Using Eq. (13.61) in Eq. (13. 4p one finds 

(5,,5,, - 5,/,j0 = (5,,((5,,0) - 5e,(5.,0) = -i{e^a^tl - t2a^e\) 9^0. (3.11) 

We have found that the commutator of two supersymmetry transformations gives us back 
the derivative of the original field. Since corresponds to the generator of spacetime 
translations P^, this equation implies the form of the supersymmetry algebra that was 
foreshadowed in Eq. (12.21) . 

For the fermion %p the commutator takes the form 

(5,2(5,, - = -i(aVi)a €28^^ + i((T^e^)a ei^^^ . (3.12) 

After applying the Fierz identity 

with X = c'^ej, C = €2, rj = df^ip, and again with x = (^^A^ ^ = V = d^^ip, followed in 
each case by an application of the identity 

e^fT'^x = -x'j^^^ = {x^(r>^o* = -(eo^V)* (3.14) 

we obtain 

{6,^6,, - 6,,6,2)ilJa = -iieia^'el - e2a^e{) di^ipa + ieia ela^df^ip - ie2a ela^'d^ip . (3.15) 

The last two terms vanish on-shell; that is, if the equation of motion a^d^ip = following 
from the action is enforced. The remaining piece is exactly the same spacetime trans- 
lation that we found for the scalar field. The fact that the supersymmetry algebra only 
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closes on-shell (when the classical equations of motion are satisfied) might be somewhat 
worrisome, since we would like the symmetry to hold even quantum mechanically. This 
can be fixed by a trick. We invent a new complex scalar field F, which does not have 
a kinetic term. Such fields are called auxiliary, and they are really just book-keeping 
devices that allow the symmetry algebra to close off-shell. The Lagrangian density for 
F and its complex conjugate is simply 



■-auxiliary 



F*F (3.16) 



which leads to the equations of motion F = F* = 0. The new field is now included in 
the supersymmetry transformation rules 

SF = -ie^a^'di.ijj , SF* = id^iP^a^'e . (3.17) 

Now the auxiliary part of the Lagrangian density transforms as 

<^>Cauxiiia.y = -is^a'^d^i^ F* + id^t/j^a'^e F (3.18) 

which vanishes on-shell, but not for arbitrary off-shell field configurations. By adding 
an extra term to the transformation law for and 

5V^„ = -i{a^e% d^ct> + e„F , 5^ = i(ea'^) . d^cp* + e^F* , (3.19) 

one obtains an additional contribution to 5>Cfermion; which just cancels with 5i2auxiiiary, up 
to a total derivative term. So our "modified" theory with C — >Cscaiar + >Cfermion + >Cauxiiiary 
is still invariant under supersymmetry transformations. Proceeding as before, one now 
obtains for each of the fields X — 0, 0*, ■0, ^0^, F, F* that the supersymmetry algebra 
now closes also off-shell 

{5J,, - = -i{e^a''el - e^a^^el) d^X . (3.20) 

The real reason for the necessity to introduce the auxiliary field F is that the two 
fields (f) and i/j have a different number of degrees of freedom. While on-shell, the two 
real propagating degrees of freedom of the complex scalar field match with the two 
spin polarization states of ip, there is a different situation off-shell. Off-shell, the Weyl 
fermion is a complex two-component object and has four real degrees of freedom (half 
of the degrees of freedom are eliminated by going on-shell). The difference of degrees 
of freedom off-shell is corrected by the introduction of two more real scalar degrees of 
freedom in the complex field F, which vanishes going on-shell. 

Invariance of the action under a symmetry transformation always implies the existence 
of a conserved current, and supersymmetry is no exception. The supercurrent is an 
anticommuting four-vector carrying a spinor index. With the Noether procedure we 
receive 

e/'. + e./t-.E^^-^" (3-21) 

where X = (f), 0*, ip, ip^ F, F* and A'^ is an object whose divergence is the variation of the 
Lagrangian density under the supersymmetry transformation, 6£ = d^j^K^. Calculated 
explicitly, and its hermitian conjugate become 

= {a^a^tl^)^ , 4^ = {iP^a^a^a . (3.22) 
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Using the equations of motion, it can be shown that the supercurrent and its hermitian 
conjugate are conserved separately: 

d,J^ = 0, 8,4^ = 0. (3.23) 

The corresponding conserved charges to these currents are 

Q^ = y/2 J d^xJ^, Qi^V2j d^xjf . (3.24) 

As quantum mechanical operators, they satisfy 

[eQ + e^Q\ X] = -iV2 5X (3.25) 

for any field X, up to terms that vanish on-shell. This further leads (using the canonical 
equal-time commutation and anticommutation relations) to the supersymmetry algebra 

{g«, Qi} = 2< . , {g„, Q^} = {Ql gt } = , (3.26) 

and since we only deal with global supersymmetry transformations we have 

[g,,P1 = o, [gt,p'^] = o. (3.27) 

So, Qa and Q\ finally can be really identified as the generators of the supersymmetry 
transformation. 

3.2 Non-Gauge Interactions of Chiral Supermulti- 

plets 

In a realistic theory like the MSSM, there are many chiral super mult iplets, with both 
gauge and non-gauge interactions. In this section, we will construct the most general 
possible theory of masses and non-gauge interactions for particles that live in chiral su- 
per multiplets. We will find that the form of the non-gauge couplings, including mass 
terms, is highly restricted by the requirement that the action is invariant under super- 
symmetry transformations. 

We start with the Lagrangian density for a collection of free chiral supermultiplets la- 
beled by an index i, which runs over all gauge and flavor degrees of freedom. Since we 
want to construct an interacting theory with supersymmetry closing off-shell, each su- 
permultiplet contains a complex scalar 0j and a left-handed Weyl fermion ■^j as physical 
degrees of freedom, plus a complex auxiliary field Fi which does not propagate. The 
results of the previous section tell us that the free part of the Lagrangian is 

Aree = O'' (f)*' d ^cj)^ + ^^'a^^^V^i + F*' F, (3.28) 

where we sum over repeated indices i (not to be confused with the suppressed spinor 
indices), with the convention that fields 4>i and tfji always carry lowered indices, while 
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their conjugates always carry raised indices. It is invariant under the supersymmetry 
transformation 

S(t)i = e^Pi, 5(f)*' = eV^', (3.29) 

5(A)a = -i{(T''e% d^(t>i + e^Fi, 5(V'^')d = i{e<j'')a 9^0" + e\F*\ (3.30) 

5Fi = -ie'^ai'd^^Pi, 5F*' = id^^^'a^'e . (3.31) 

The most general renormalizable interaction Lagrangian of these fields that is invariant 
under supersymmetry transformation can be written as 



'int 



(^-hw'j^.^. + w'F^j + c.c. (3.32) 



with the so called superpotential W (which is not a potential in the usual sense) 

W = ^M^,(P, + ^y''''Mj<Pk , (3.33) 

6W ■■ 6^ 

d(pi d(pid(pj 

M*^ is a symmetric mass matrix for the fermion fields, and y"^^^ is a Yukawa coupling 
of a scalar 0^ and two fermions ipiipj that must be totally symmetric under interchange 
of i,j,k. So we have found that the most general non-gauge interactions for chiral 
supermultiplets are determined by a single analytic function of the complex scalar fields, 
the superpotential W. 

The auxiliary fields Fi and F** can be eliminated using their classical equations of motion. 
The part of ^frcc + 'Cint that contains the auxihary fields is FiF*' + Fi + W* F*\ leading 
to the equations of motion 

Fi = -W* , F*' = -W' , (3.35) 

so the auxiliary fields can be expressed by terms corresponding to the superpotential. 
After making this replacement, the Lagrangian finally takes the form 

Ahirai = d'^^d^cPi + iiP^'a^d^^Pi - i (PF^ViV^,- + W*^,pH^^) - W'W* . (3.36) 

Another way of writing the Lagrangian is to divide it into kinetic-, mass- and Yukawa- 
couphng terms and a scalar potential y(0, 0*) 

>Cchirai = d^<t>*'d^4>i - V{(t>, (t>*) + i^p^'a^d^^Pi - ^M^'i^ii^j - ^M^-V^t^^b' 

—y'^'Mji^k - \y*iur^^^^^^ (3.37) 



with 



1/(0, 0*) = WWl = F*''Fk = 

M;^M^^rh + \M"^y%^,(f'<f' + \M:^^"r'h<t>k + i|/^^"i/k0.0.0*V*' • (3.38) 
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Now we can compare the masses of the fermions and scalars by looking at the hnearized 
equations of motion: 

d^d^<f>i = -M*^M''^<j>^ + ..., d^d^A = -M^M'^V,- + . . . . (3.39) 

Therefore, the fermions and the bosons satisfy the same wave equation with exactly the 
same squared- mass matrix with real non-negative eigenvalues, namely (M^)/ = M*j.M''^ . 
Since SUSY is not broken yet, we have a collection of chiral supermultiplets, each of 
which contains a mass-degenerate complex scalar and Weyl fermion. 

3.3 Lagrangians for Gauge Supermultiplets 

Besides the chiral supermultiplets there also exist gauge supermultiplets that contain 
gauge bosons and their fermionic superpartners. The gauge bosons are described (be- 
fore spontaneous symmetry breaking) by a massless gauge boson field A^, the superpart- 
ners by a two-component Weyl fermion gaugino A". The index a runs over the adjoint 
representation of the gauge group (a = 1, ... ,8 for SU{3)c color gluons and gluinos; 
a = 1,2,3 for SU{2)l weak isospin; a = 1 for U{1)y weak hypercharge) . The gauge 
transformations of the vector supermultiplet fields are 

S„A; = -a^A" + gr'^A^A^ , (3.40) 
W^" = (3.41) 

where A" is an infinitesimal gauge transformation parameter, g is the gauge coupling, 
and f"''"^ are the totally antisymmetric structure constants that define the gauge group. 
The on-shell degrees of freedom for A^ and amount to two bosonic and two fermionic 
helicity states (for each a), as required by supersymmetry. However, off-shell A^ con- 
sists of two complex, or four real, fermionic degrees of freedom, while only has three 
real bosonic degrees of freedom. So, we will need one real bosonic auxiliary field, tra- 
ditionally called D", in order for supersymmetry to be consistent off-shell. This field 
also transforms as an adjoint of the gauge group and satisfies {D"-)* — D"-. Like the 
chiral auxiliary fields Fj, the gauge auxiliary field has no kinetic term, so it can be 
eliminated on-shell using its algebraic equation of motion. 
The Lagrangian density for a gauge supermultiplet is thus given by 

Aauge = ^\f%F^'''' + i\^^a^D,\^ + ^L>«L>" , (3.42) 

where 

= d,Al - d^Al - gr^'^A^Al (3.43) 
is the usual Yang-Mills field strength, and 

D^X^ = 8^X^-9 r^'^A^X^ (3.44) 

is the covariant derivative of the gaugino field. Up to multiplicative factors, the super- 
symmetry transformations of the fields are 

M« = -^(eta^A'^ + At%6) , (3.45) 
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= -^(^'^'^^)" + 7f ^'^ ' (3.46) 

= — ]= {t^af'D^X - D^X^'^a^'e) . (3.47) 
v2 

Like the transformations of the chiral fields, they satisfy the commutator relation 

{6,A. - KSe,)X = -t{e,a^el - e2a^e\)D^X (3.48) 



for X equal to any of the gauge-covariant fields F^^, A*^, A^", D"^, as well as for arbitrary 
covariant derivatives acting on them. This ensures that the supersymmetry algebra is 
realized on gauge-invariant combinations of fields in gauge supermultiplets, as they were 
on the chiral supermultiplets. 

3.4 Supersymmetric Gauge Interactions 

Finally we are ready to consider a general Lagrangian density for a supersymmetric the- 
ory with both chiral and gauge supermultiplets. Suppose that the chiral supermultiplets 
transform under the gauge group in a representation with hermitian matrices {T°')i sat- 
isfying [T",T''] = if'^^T". [For example, if the gauge group is SU{2), then = e"-^", 
and the T'* are 1/2 times the Pauli matrices for a chiral supermultiplet transforming 
in the fundamental representation.] Since supersymmetry and gauge transformations 
commute, the scalar, fermion, and auxiliary fields must be in the same representation 
of the gauge group, so 

= W^^'iT'^X), (3.49) 

for Xi = (pijipi, Fi. To have a gauge-invariant Lagrangian, we now need to replace the 
ordinary derivatives in Eq. fl3.28p with covariant derivatives: 

d^<P, -> D^<l), = d^(l), + tgA'^^{Ty), (3.50) 
^ D^0« = 9^0«-^(7A;:(0*n* (3.51) 
d^^P, ^ D^^lji = d^^, + igA';^iT''^)i. (3.52) 

This simple procedure achieves the goal of coupling the vector bosons in the gauge super- 
multiplet to the scalars and fermions in the chiral supermultiplets. Nevertheless, there 
are three additional possible interactions that are gauge invariant and renormalizable 
and can be included in the Lagrangian, namely 

(0*T»A^ A^'^(^/'^T», and (0*T»D^ (3.53) 

It is only possible to add these terms to the Lagrangians for the chiral and gauge super- 
multiplets, if the supersymmetry transformation laws for the matter fields are modified 
to include gauge-covariant rather than ordinary derivatives. Also, it is necessary to 



include one strategically chosen extra term in 6Fi, so 

S^i = etPi (3.54) 

Sipia = -i{a^'e^)aD^(f)i + eaFi (3.55) 

6Fi = -ze^a'^D^^, + V2g{Ty),e^X^\ (3.56) 
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Like the auxiliary fields Fi and F**, the are expressible in terms of the scalar fields, 
using an equation of motion 

D"" = -g{<i)*Ty) . (3.57) 
One thus finds that the complete scalar potential (see Eq. fl3.38p ) becomes 

1^(0, 0*) = + i 5^ D'^D'^ = W:W' 9l{<P*Tyf . (3.58) 

a a 

Summing up all the previous results we obtain as the full Lagrangian density for a 
renormalizable supersymmetric theory 

•^SUSY -^chiral ~l~ -^gauge 

-V2^(0*T»A" - V2gX^%^^T^<f)) + g{(j)*T^(j))D'' (3.59) 

-\Plg ((0*T»A'^ + A^''(^/'^T») . (3.60) 

The first line of Eq. fl3.60p shows the kinetic terms of fermions and scalars, the self- 
interaction of gauge-fields and the kinetic term of gauginos. In the second line, the first 
term is the scalar-potential while the remaining terms come from the superpotential 
and include Yukawa-couplings and fermion-mass terms. The last line consists of the 
additional supersymmetric couplings. 

3.5 Soft Super symmetry Breaking Interactions 

Since no supersymmetric particles are discovered yet, they must have significantly higher 
masses than the standard model particles. In other words, supersymmetry is broken. 
Like electroweak sjTiimetry breaking, supersymmetry should be broken spontaneously. 
Although this mechanism of producing the symmetry breaking at high energies is not 
totally understood yet — and there is no consensus on which of the several models is the 
right choice — it is possible to write down the general form of possible supersymmetry 
breaking terms in the Lagrangian at low energies 

/:soft = - Qm, A'^A" + ia^^ V*0,0fc + + c.c. - (m2);.0^*0, (3.61) 

with gaugino masses Ma for each gauge group, scalar squared- mass terms ivp?y^ and 
and (scalar)^ couplings a*-'^. 

Because the effective Lagrangian can be written in the form L = Csvsy + -^soft (see 
the next section) and because £soft contains only mass terms and coupling parameters 
with positive mass dimension, broken supersymmetry is still providing a solution to the 
hierarchy problem. It is therefore named "soft" supersymmetry breaking. 
Furthermore, supersymmetry is indeed broken by £soft; because it involves only scalars 
and gauginos and not their respective superpartners. 
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3.6 The Complete Lagrangian 

Finally, after all the efforts made in this chapter, we can present the full Lagrangian of a 
super symmetric theory, consisting of the full supersymmetric Lagrangian from Eq. (13.601) 
and the soft supersymmetry breaking Lagrangian from Eq. (13.611) resulting in 

1 



1 
6 



\Ma A" A" + \a'^^(t)i(t)j(t)k + \h''(t)4i ) + c.c. - (m^);.^^*^^ . (3.62) 
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Chapter 4 



The Minimal Supersymmetric 
Standard Model (MSSM) 

The MSSM extends the particle content of the Standard Model in two different ways [HI 
[20] . On the one hand, all particles get a superpartner, on the other hand there is a 
larger Higgs sector with two complex Higgs doublets. The superpartners get the same 
names as their corresponding SM-particles, just with the prefix "s" (scalar) for spin = 
superpartners and the suffix "-ino" for spin = 1/2 superpartners. That leads to names 
like sfermions on the one, and gauginos and higgsinos on the other hand. Particles and 
superpartners are placed in the chiral and gauge supermultiplets as shown in Table 14.11 
and 1121 



Names 


spin 


spin 1/2 


SU{3)c, SU{2)l, U{1)y 


squarks, quarks 


Q 


{ul di) 


{ul di) 


(3, 2, I) 


(x3 families) 


u 






(3, 1, -I) 




d 


d* 


d^R 


(3, 1, i) 


sleptons, leptons 


L 


('^ ei) 


ei) 


(1, 2, -i) 


(x3 families) 


e 


eR 


e^R 


(1, 1, 1) 


Higgs, higgsinos 


Hu 


{Hi Hi) 


m Hi) 


(1, 2, +i) 






{Hi Hi) 


{Hi Hi) 


(1, 2, -i) 



Table 4.1: Chiral supermultiplets in the MSSM. The spin = fields are complex scalars 
and the spin = 1/2 fields are left-handed two-component Weyl fermions. 

The left-handed and right-handed pieces of the quarks and leptons are separate two- 
component Weyl fermions with different gauge transformation properties in the SM, so 
each must have its own complex scalar partner. It is important to keep in mind that 
the "handedness" here does not refer to the helicity of the sfermions (they are spin-0 
particles) but to that of their superpartners. 
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Names 


spin 1/2 


spin 1 


5f/(3)c, SU{2)l, U{1)y 


gluino, gluon 


9 


9 


( 8, 1 , 0) 


winos, W bosons 


A± A3 




(1, 3, 0) 


bino, B boson 


A' 




(1, 1, 0) 



Table 4.2: Gauge supermultiplets in the MSSM. 

The symbols for the superpartners carry a tilde (~) for distinction. 
Due to the spontaneous electroweak symmetry breaking, the interaction eigenstates 
(winos, bino) are no longer mass eigenstates and thus no physical particles. Therefore 
one needs to consider the mixing of interaction eigenstates to mass eigenstates as we 
will do in detail in Chapter [51 Table 14.31 shows the mass eigenstates and corresponding 
interaction eigenstates. 



Names 


Spin 


Pr 


Gauge Eigenstates 


Mass Eigenstates 


Higgs bosons 





+1 


Hi Hi Ht H^ 


H^ A^ H^ 


squarks 





-1 


ul ur di dn 
cl cr Sl Sr 
ti tR bi bR 


Ul U2 di d2 
Cl C2 Sl S2 

ti t2 bi 62 


sleptons 





-1 


Cl Cr Ve 
Tl Tr Vr 


ei 62 Ve 
/Ul /I2 z/^ 
Ti T2 Vr 


neutralinos 


1/2 


-1 


A' A3 Hi Hi 


Xl X2 X3 X4 


charginos 


1/2 


-1 


A± Hi Hi 


x^ x^ 


gluino 


1/2 


-1 


9 


9 



Table 4.3: Mass eigenstates and corresponding interaction eigenstates of the particles in 
the MSSM. 



Chapter 5 
Mass Matrices 



Due to the spontaneous electroweak symmetry breaking, the mass eigenstates (and there- 
fore the physical particles) of the superpartners are not identical with the gauge eigen- 
states of the interaction, but mixtures of them. Thus one has to deal with the mixing 
of the fields and derive mass matrices, which can then be diagonalized using rotation 
matrices in order to obtain mass eigenstates. Since these rotation matrices include in 
general complex values, they are an important source of CP violation. 



5.1 Sfermion Sector 

The relevant terms for the sfermion mass matrix are derived from the soft SUSY breaking 
terms and the auxiliary field terms (for a detailed derivation see |15]). For the mass 
matrix in the basis ip = (/l /j?)''' (/ = {t, b,f, . . .}) we get 

Mi = ( '"^^ ] = I (5.1) 

I ™Ri '"rk / I '"/«/ in' 



with the following entries (see also [2U]): 



mj-^ = Mfg.^y + mj + m| cos 2/?(j|^ - e/ sin^ Ow) , (5.2) 

,2 ]\4'2 



my^ = M^-.^.-y + rrif + cos 2(3 ej sin 6w , (5.3) 

r m„(A* — yU cot /9) ... up-type sfermions 
m/O/ = < (5.4) 
^ rriii^A*^ — yUtan/3) . . . down-type sfermions . 

M^Q-Lj, Mjjj.j.g}. are real soft SUSY breaking masses, A^, are complex trilinear break- 
ing parameters and is the complex higgsino mass. Ij^ is the third component of the 
weak isospin, e/ is the electric charge in terms of the elementary charge, rrif, rriu, rrtd are 
the masses of the fermionic superpartners and tan /? = ^ is the ratio of the two VEVs 
of the Higgs fields. 

The mass matrix can now be diagonalized with a rotation matrix 

I w?- nifaf \ r . I rri'] \ ? 



I / \ m| 
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R{r ^ 


H 


[ rL 


rL ) 





We parameterize the unitary matrix as 

f I -"-ir -"-iR \ I cos6'r- e*'^/sin6'r . 

-e sm Uj cos U j 

with the mixing angle 6j and the complex phase (pj which is a new source for CP 
violation. 

Applying the rotation matrix to the gauge eigenstates /„ (a = L, R) one obtains the 
mass eigenstates fi{i = 1,2) and vice versa: 

_ roAt/^ /A ^ f _ of* 



IrJ ^^^^ V/J ^ f- = ^f^^ (5.8) 

with the convention m^^ < mj^. We assume that the super CKM matrix is diagonal, so 
there is no mixing between the three generations of sfermions. 



5.2 Chargino Sector 

The charged fermionic superpartners and , ip'jj^ of the W and Higgs bosons mix 
to two mass eigenstates with charge ±1, the so called charginos. We identify these 
eigenstates in the Dirac spinor representation with xf {i = 1, 2) using the convention 
m-± < m-± (see [H], [16], [20]). 

The relevant terms for the chargino mass matrix are derived from the soft SUSY breaking 
terms, the SUSY gauge interaction terms and the Yukawa interaction terms, which can 
be deduced from the superpotential (for a detailed derivation see [45]). 
The mass term of the lagrangian in the basis = {—iX'^ "^//a)'''' ^~ ~ {—iX~ 'ipn-^)^ is 
given by 

= -larVirV) -(^l ""^y (^') + h.c. (5.9) 
with the mass matrix 

v2mvycosp iJ, 

M2 is the real SUSY breaking mass of the gauginos and A^, the superpartners of the 
SU{2)l bosons A^ and A"^. The parameter fi is the complex higgsino mass. 
This matrix can be diagonalized with two 2x2 rotation matrices U and V: 

, m~± 

Xd = U*XV-' = I ^1 I . (5.11; 

"^xt 
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Applying the unitary matrices U and V to the gauge eigenstates ^pf one obtains the 
mass eigenstates xf and vice versa: 



XT 



^7 = U*,x. 



Finally we can define the charginos as Dirac spinors {i = 1,2): 



(5.12) 
(5.13) 



'X.i 



'X.i 

xt 



(5.14) 



5.3 Neutralino Sector 

The uncharged fermionic superpartners of the f/(l)y interaction field B^, the third 
component of the SU{2)l interaction field A^^ and the Higgs bosons Hi and if| mix to 
four mass eigenstates called neutralinos. We specify them in the Majorana representation 
with x° = 1, 2, 3, 4) using the convention m^^o < m^o < m^o < m^o (see [H], [16], [20]). 



The relevant terms for the neutralino mass matrix are again derived from the soft SUSY 
breaking terms, the SUSY gauge interaction terms and the Yukawa interaction terms, 
which can be deduced from the superpotential (for a detailed derivation see [IS]). 
The mass term of the lagrangian in the basis ^/'^ = {—iX', —iX^, V'lfa) given by 



0\T 



with the mass matrix 
/ 



Y 



Ml 




\ mzse^sp 





M2 



+ h.c. 



-mzse^cp 
mzce^cp 




mzse^sp 
-mzce^.sp 





(5.15) 



(5.16) 



with the abbreviations '■= sin a and := cos a. 

The matrix can be diagonalized with the unitary rotation matrix Z: 



Y, 



D 



Z*YZ- 



diag(m^o,m^o,m^o,m^o) . 



X3' 



X4 ^ 



(5.17) 



Applying the matrix Z to the gauge eigenstates ip^ one receives the mass eigenstates 
and vice versa: 

X° = ^.,V'° V^° = ^*X°. (5.18) 
At last the neutralinos are defined as Majorana spinors with mass index {i = 1, 2, 3, 4): 
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Chapter 6 



Couplings 



In this chapter we provide the relevant terms of the interaction Lagrangian (and thus 
the couphngs) for the leading contributions mentioned in Chapter [9] and some more. For 
a listing of more couplings see [201 US HZ] • 

Note that although we have written the fermions and sfermions only in the third gen- 
eration we sum over all generations nevertheless. We apply the Einstein summation 
convention, repeated indices are summed. We use the same convention as in pO] . 
The chargino-squark- quark {xtl^'i) chargino-slepton-lepton (xt^^'i) interaction 

Lagrangian can be written as 



+xt{Aj::PR + A^*PL)th* + xp{Bt*PR + B^:PL)bi* (e.i) 

= MAi'^PR + <PL)xtr^ + fiBi^PR + Bi^PL)xru, 

+n{A'tPR + A',fPL)^rf: + WiBtPR + B'tP,)Ti>: (6.2) 



with the projection operators Pr,l = (1 ±7^)/2. The abbreviated coupling matrices are 




(6.3) 



(6.4) 
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The neutralino-squark-quark (XkQQi) the neutrahno-slepton-lepton {x^lh) interac- 
tion Lagrangian can be written as 

>CxO,,- = ^C^,PH + C^iPL)& + b{D^n + D^,Pr.)xlk 

+xUCt*PR + C^*PL)ti* + xI{DI:Pr + D^:P,)hh* (6.5) 

+WC'u'*Pr + C',^*PL>rK + Md'^:Pr + ^r^L)rf* . (6.6) 
The abbreviated couphng matrices are 



Dki = —?^\ o^k^Rii + {Zk2 - -iiiIl9wZkl)Ril 

C'lf ~ (— '^A;2 + tan6'vK-^iki) 
= 



= 7! - . (6.8) 

The gluino-squark-quark {gqqi) interaction Lagrangian is 

= WfPR + G^PL)qq* + ?(G'f *Pii + Cf ^P^)^ (6.9) 
with the couphng matrices 

Gf = V2^.TP| 

Gf = -V2g,TRl,. (6.10) 

gfs is the couphng constant of the strong interaction and T is the generator of the SU {3)c 
group. We omitted the colour indices of the quark and squark and the gluino index for 
the sake of simphcity. 

The photon-chargino-chargino {jxlxi) interaction Lagrangian is 

£,;e+x+ = eAf'^^^xt (6-11) 
with the electric positron charge e. 

The Lagrangian of the photon-sfermion-sfermion {'jfifi) interaction can be written as 

C^jj=-ieefA>^f:d,f, (6.12) 

where e/ is the charge of the particle in e (e.g. for the electron e/ = — 1). If we define 
the momenta fi{ki) and //(/C2) we obtain as the Feynman rule — iee/(/ci — A;2)^- 



Chapter 7 
Definitions 



In this chapter we give some definitions which are important for the later study of 
CP violating effects. After a brief information about the CP transformation (see [15] ) 
we define the CP violating decay rate asymmetry 6'"^ (see [20], for more details see 
Chapter [8]). Further we define the branching ratio BR and explain why it is important 
to observe both 6'"^ and BR simultaneously. 



7.1 CP Transformation 

The C transformation (charge transformation) changes particles to its anti-particles and 
vice versa. It changes the sign of the electric charge and inverts the colour charge of a 
particle. For spinors we have the relations 

u{k) = Cv'^ik) u{k) = v'^{k)C 

v{k) = Cu^{k) v{k) = u~^{k)C (7.1) 

and for the vector bosons 

< = Ce, . (7.2) 

The P transformation (parity transformation) is equivalent to a spatial point refiection. 
It changes the sign of the coordinate system and thus a right-handed coordinate system 
into a left-handed one or vice versa. For spinors we use the definition 

Pip{t,x)=-f^ip{t,-x). (7.3) 

Together they form the CP transformation. Applied to the Lagrangian, it changes 
signs of four momenta, left- and right-handed parts and spinors and therefor it also 
conjugates the tree-level coupling matrices. Since these are now complex (due to the 
complex chargino rotation matrices U and V, the complex neutralino rotation matrix 
Z and the complex sfermion rotation matrix R) we have a violation of CP symmetry. 
As we assume the Cabibbo-Kobayashi-Maskawa (CKM) matrix (and the super CKM 
matrix as well) to be diagonal we neglect the small CP violating effects coming from 
flavour mixing. 
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7.2 Decay Rate Asymmetry S^^ 

The CP violating decay rate asymmetry 6'"^ for our decay is defined by 

with the decay widths 

r+ (X $^IM+eel' + 2Re(5^(M+J^-MLp) (7.5) 

s s 

T oc 5^|M-,,j2 + 2Re(5^(M-jt>lr„„p). (7.6) 

s s 

The matrix elements at tree- and one-loop-level are given by 

7W+ee = tu{kl){B''PR + B''PL)v{~h) 
-Mt'ree = t uih)iB''* Pr + B^^* PlM-^) 

= tuih)i5B^Pn + SB'^PLM-h) 
M^op = tu{k2){SB^Pn + 6B'lPL)v{-h). (7.7) 

Since there is no CP violation at tree- level, we set |A^treeP := l-^treeP = l-^treeP- The 
decay rate asymmetry can then be approximated to 



9r 

tree 

2Re(E.(^iee)^A^ltop) - 2Re(E.(^tree)^A^roop) 



loop J 

with 



2EJMree| 



(7. 



2Re(5^(MU^MU) = 2{ml-ml-ml,)Re{B^*6B^ + B^*6B^^) 

s 

-Ambm^+Re{B^*6Bl + B^*5B^) (7.9) 
2Re(5^(M-jt>lr„„p) = 2(m|-m^-m^,)Re(i?W + i?^<55^) 



-Ambm-+Re{B^6B^ + B^5B^) (7.10) 

and 

2 5^ lA^treeP = 2(m| - - m2^)(|i?^|2 + _ 8m,m-+Re(S^*i?^) . (7.11) 

s 

In Chapter we will only provide the form factors of the graph but not the anti-graph, 
so Aji^L = {4:iiySB^'^ . The form factors 6B^'^ of the anti-graph can be easily obtained 
by conjugating all the couplings involved. 
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We can refine our result by deriving more specific expressions for Eq. (17.91) and (I7.10p . 
We rewrite these equations to 

<5r±cx2Re(^(A1it,,Jt^±„p) = 2ARe{B^6B2 + B^6B^) 

s 

-Arrn,m^+Re{B^5B^ + B^5Bl) (7.12) 

using A = (m? — — ""^^h-), -B+'^ = B^'^ and B^'^ = B^'^*. Then we define the 
combined coupling matrices 

C'l = Bl6Bi Cl = B\5Bi (7.13) 
with i,i e {i?, L} and we get 

5T^ oc 2 A (Re(Cf^) + Re(C|^)) - 4m6m-+ (Re(Cf^) + Re(C|«)) . (7.14) 

The coupling matrix C^? can be generally expressed by (see Eq. (19.21) and (19. 3p for 
specific examples of 5B±) 

C% oc 6 X ^o^i^2 X PaVe oc b* x (gogm)* x PaVe (7.15) 

where b = b* = 5+ is the coupling at tree-level, gogig2 are the couplings of the three 
vertices and PaVe are the Passarino-Veltman-Integrals (see Chapter [9] and Appendix iBl). 
Taking the real part of the coupling matrix Cj? and keeping in mind that both the 
couplings and the PaVe's are in general complex we get 

Re(6^o^i^2 X PaVe) = Re(6^o^i^2)Re(PaVe) - Im(6^o^i^2)Ini(PaVe) (7.16) 
Rei{bgogig2y x PaVe) = Re(6(7o^7i^?2)Re(PaVe) + Im(6^o^i^2)Im(PaVe) . (7.17) 

This finally leads to the important decomposition into both CP invariant and CP vio- 
lating parts 

^iCl) = Ct±lcg, (7.18) 

with the definitions 

Ct oc Re(6^o^7i^2)Re(PaVe) (7.19) 
Ci^P oc -2Im(6(7o^7i^72)Im(PaVe) . (7.20) 

To show that the decomposition in Eq. (I7.18P is correct one simply has to calculate 
S^^ oc 6T^ — 6T~ using the rewritten equation 

6T^ oc 2A(Qt^ + C^,^±i(C^« + 0) 

-4m,m,. (c^J; + Ct^ ± \{C§^ + CS^)^ (7.21) 

resulting in 

6^^ oc 2A(C^^ + C^^) - AmhTn^+ {C^^ + C^^) . (7.22) 

From Eq. (I7.20p we can see that we need not only the couplings but also the PaVe's to 
be complex in order to obtain a non-zero S'"^ . 
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7.3 Branching Ratio BR 

The branching ratio BR is simply defined as the ratio between the tree-level decay width 
of a certain decay and the total tree-level decay width. We write 

BR-^r- (7.23) 

tree 

An example would be 

BR{U - bxt) = "^0^ . (7.24) 

ii 

This quantity gives you the probability how often the certain decay occurs compared 
to all possible decay channels. Combined with the decay rate asymmetry to 6"^^ x BR 
one obtains the overall probabihty how often the certain decay channel is CP violated 
compared to the rest. As S*"^ usually gets high when BR gets low and vice versa one 
has to find an optimum in S"-^^ x BR for a good measurability at colliders. 



Chapter 8 
CP Violation 



In this chapter we provide a deeper insight into CP violation by stating more precisely 
how and when the CP violating decay rate asymmetry 6^^ defined in Chapter [7| occurs. 

As we already found out, the tree- level couplings introduced in Chapter [6] being 
complex is just a necessary but not sufficient condition for CP violation. There are in 
total two conditions which have to be fulfilled in order to obtain a non-zero decay rate 
asymmetry 6'"^ at one-loop level: 

1. Complex tree- level couplings 

2. (a) At least two open decay channels 

or equivalent ly 

(b) At least one Passarino-Veltman-lntegral has to become complex 
or equivalent ly 

(c) At least two particles in the loop have to become on-shell. 

The first condition can be easily understood by remembering that a CP transformation 
of a Lagrangian results in the conjugation of the tree-level couplings involved (see Chap- 
ter [7]). Real tree-level couplings leave the CP symmetry intact. Pragmatically one can 
also see this condition in Eq. f l7.20p . 

Furthermore, the reason why 6'"^ does not occur at tree-level can be seen by looking at 
Eq. (17.111) . keeping in mind that F oc |A^p = '^^A4^A4. The tree-level couplings 
B^'^ and their complex conjugates B^'^* appear symmetric in the equation, so a CP 
transformation (i.e. conjugating the couplings) keeps the decay rate CP invariant. An 
other way to see this is the fact that both CP transformation and calculating an adjoint 
matrix conjugate the couplings so the net effect is zero. The decay rate asymmetry 
is thus a pure loop-effect starting at one-loop level. 

For the second condition we need to expatiate on the subject in order to understand the 
three conditions and their equivalency. We start with the requirement 2. (a). Following 
from the CPT Theorem one can prove that the total decay width Ftotai of a certain 
particle is invariant under CP transformation. That means if we we assume at least two 
different kinematically possible decays we have 

rt+otai = r+ + F+ = F- + F- = F-,,, . (8.1) 
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Note that only the total decay width keeps CP invariant but the partial decay widths are 
not necessarily CP invariant! So in our example we have in general P^ 7^ Fj", P^ 7^ P^ 
and we can write 6'-^^ oc — ^12 0- Only in the case of just one open decay channel 
the partial decay rate (being now the total decay rate) is becoming always CP invariant. 

The requirement 2.(b) can be directly seen in Eq. fl7.20p . If Im(PaVe) = the CP 
violating part becomes zero thus rendering 6'"^ in Eq. f l7.22p zero as well. 
The last requirement 2.(c) can be understood in two ways. First one can calculate the 
discontinuity (and hence the imaginary part) of a PaVe using Cutkosky Rules ^911501 151]. 
After applying these rules one directly obtains the condition that the particles in the 
loop have to become on-shell for a non-zero Im(PaVe). For example, if we apply the 
rules and cut through the loop of the Generic Structure I in Chapter [9], we receive the 
condition rriQ > Mi + M2 for Im(PaVe) 7^ 0. If we take the Generic Structure II we 
can cut in three ways (vertically and two times horizontally) and get the conditions 
mo > Ml + M2, mi > Mo + Mi or m2 > Mq + M2 for Im(PaVe) ^ 0. In both examples 
we need at least two loop-particles who have to be on-shell, which is exactly the require- 
ment 2.(c). Note that the condition mj > Ylj^j implies that the decay of a particle i 
into other particles j is kinematically possible, i.e. r{i j) 7^ 0. 

The second way to understand the last requirement 2.(c) is by using the Optical The- 
orem [51j directly, which we will derive very briefly. From unitarity of the S-matrix 
'^kSji.Ski = Sfi and the relation Sfi = 6fi + iAiji to the matrix element A4 we can 
write 

k 

= 6fi-iM}, + iMf, + J2i^}kMki) (8.2) 

k 



with the indices i, /, k denoting initial-, final- and intermediate-states, respectively. 
Rewriting this condition directly results in the Optical Theorem 

2lm{M /,) = J2(^}k^k^) oc atotal (8.3) 
k 

which relates the scattering amplitude A4 fi to the total cross section (Xtotai of the process. 
Note that if we set i = f the theorem relates J^a with the total decay width Ptotai oc 
J2k l-^fejP of the process i —>■ k. Without loss of generality, we choose the i k process 
to be a tree- level decay Ptree(^ k) of one particle i into two others. The amplitude J^a 
then represents a one-loop selfenergy of the particle i. Now the theorem states that the 
imaginary part of the scattering amplitude M-u (and hence the imaginary part of the 
PaVe inside of M-u) becomes non-zero only in the case of a kinematically possible decay 
rtree(^ k) , i.e. iTii > Ylik^k- Pictorially speaking, the one-loop selfenergy of particle 
i can then be "assembled" by two tree-level processes (or equivalently "cut" through to 
get two tree-level processes) as shown below: 
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This is exactly what Cutkosky Rules are doing and so these two ways to understand the 
requirement 2.(c) are equivalent. 

Because of the connection between Im(PaVe) and rrii > Mj the requirement 2.(c) is 
equivalent to requirement 2.(b). To show the equivalency with the first requirement 2. (a) 
we take the Graph 1 in Chapter [9]as example. Let's assume that only two decay channels 
are open: U hxl and U tg. Condition 2. (a) as well as the conditions 2.(b-c) state 
that CP violation is possible. Now let's assume that only one decay is kinematically 
allowed: ti — > hxl- Condition 2. (a) states that CP is not violated. Because of mi. < 
rrit + rrig (and hence Im(PaVe) = 0) the conditions 2.(b-c) predict no CP violation as 
well. So we demonstrated that all three requirements (a-c) of condition 2. are really 
actually equivalent. 

Concluding we want to mention that selfenergy-loops like the one in Graph 1 can only 
contribute to 6'--^^ in the case of different particles, i.e. i ^ j. If the particles are the 
same, the couplings on the two vertices of the loop are complex conjugates of each other 
{9i = 92*y 9i — 92*) resulting in an CP invariance of the form factors of the loop. 
This is also the reason why there are only three different topologies in Chapter [9] who 
can contribute. The possible fourth topology, a selfenergy-loop for the particle with 
momentum ki, cannot contribute, because this particle is associated with a bottom- 
quark. As we have set the CKM matrix to be diagonal, the b-quark cannot change into 
another quark of the first or second generation, thus it has to stay the same and cannot 
fulfill the requirement. 
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Chapter 9 
Contributions 



In this chapter we hst the most important processes at one-loop level who are expected 
to yield the highest CP violating asymmetries. We start by specifying all possible 
topologies and the convention of the four momenta and the masses. Then we insert 
scalar and fermionic fields in order to obtain the most important generic contributions 
and give the matrix element and the corresponding form factors. In the last step we list 
the most important processes with their specific form factors and mention the remaining 
contributions. 



9.1 Topologies 

There are in total three possible topologies who can contribute to a CP violating asym- 
metry: two different topologies with a selfenergy contribution and one vertex correction. 
The masses of the external (mo, mi, 1712) and internal particles (Mq, Mi, M2) and the 
four momenta of the external particles [p, ki, ^2) are defined as shown below. The 
on-shell relations = ml, kf = ml, k\ = m\ as well as the relations p = ki — k2 and 
k = ki + k2 are used. 

For the calculation of the loop integrals we use the formalism of Passarino-Veltman- 
Integrals defined in the Appendix [Bl The convention of the arguments of the B- and 
C-functions are shown below. All the above conventions for momenta and masses are 
given in ^5] . 
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mo 



P 



q + k 




q + k 



mi 



m2 



= C{mlmlmlMlMlMi) 



mo 



P 




= B{mlMlMl 



9.2 Most Important Generic Structures 

Now we insert scalar and fermionic fields into the topologies in order to obtain the 
leading generic contributions. We list the matrix element and the corresponding form 
factors which can then be used for the calculation of specific processes. All form factors 
are calculated in the SUSY invariant Dimensional Reduction regularization scheme (DR) 
(see [521153]). 

The general matrix element can be written as (see [i6] ) 



M 



(477) 



(9.1) 



We only give the form factor A/j since A^ can be easily obtained by exchanging right- 
and left-handed couplings {g^ g^). The argument set of the form factors is always 
/ = (mo, mi, m2, Mq, Mi, Ms, g^, g^, g{, g^, g^). 
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9.2.1 Generic Structure I 



The first leading generic structure is a scalar selfenergy-loop with a fermion-fermion pair 
shown below. 



00, "^0 



Ml 

Mo 










lig^Pn + gl^PL) 
t{g^Pn + gfPL) 
i{g^PR + gJ^PL) 



The appropriate form factor is 

1 



ml - Mi 



9^ 



2M,M2[g^g^ + g^g^)Bo 
+ i9i92 + 9i92) (mm',) + AoiM^) + (M^ + M^ - ml)Bo 



■ (9.2) 



9.2.2 Generic Structure II 

The second leading generic structure is a scalar-fermion-fermion vertex correction. 

fi,mi 



5o,"^o 








tigl^PR + gl^Pi) 
t{g^PR + gfPL) 
^{9§Pr + 9^Pl) 



f2,m2 



Here the respective form factor is 



[9o92i9im, + g^M,) + g^g^^gi^m^ + g^M2))rmC, 
+ [9o92i9imi + g^Mi) + g^gf^g^m^ + g^M2))m2C2 
+gH 9192^1^2 + gfg^MiM^ + g^g^niiM^ + g^g^m2Mi )Co 



+^oV^72'' ( 5oK, Mi^, M|) + M^Co 



(9.3) 
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9.3 Most Important Processes 

Finally we list the most important processes with their specific form factors. They 
are the two possible processes involving a gluino g. Because the gluino couples like its 
superpartner the gluon with the strong interaction force, these contributions are expected 
to dominate over all others, if the decay channel ti ^ tg is open (m^. > rrit + nig). For 
a listing of all contributions at full one-loop level (including these leading contributions) 
which we calculated see Appendix |X1 The coupling matrices we used are defined in 
Chapter O 

9.3.1 Graph 1 

The first leading process is a stop selfenergy with a gluino-top in the loop. Note that 
i 7^ j in order to contribute to CP violation. 




The form factor has the following arguments: 

= AUK>"^^'^x^^t>^s,^^,^S,5fe,G'f^Gf^Gf,G^ . (9.4) 

9.3.2 Graph 2 

The second leading process is a vertex correction with a sbottom-gluino-top in the loop. 
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The form factor takes the following set of arguments: 

2 

^%,L = EA^^^K,m,,m-+,mg^.,m5,m,,Gf^Gf ,Gf ,Gf , Ag, . (9.5) 
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Chapter 10 

Numerical Results and Conclusions 



In this chapter we finally study the decay rate asymmetry 5^^ defined in Eq. f l7.4p and 
the combined quantity 5^^ x BR (branching ratio BR defined in Eq. (17.241) ) of the 
processes ti hxt and t2 ^ bxt (and their CP transformed counterparts) by varying 
the parameters Mq, tan/5, ipAt: ^2, A* and v?^- 

10.1 A Typical Scenario 

We fix the input parameters of the MSSM and the SM by the choice given in Table [T0?T1 
The coupling of the strong interaction force is taken running in the dimensional 
reduction regularization scheme DR [521 ESI EH ESj, renormalized at the scale of the 
decaying particle mass m^. {i = {1, 2}) in the SPA convention [56]. The gluino mass rrig 
calculated from ag via GUT relations is therefore running as well (calculated iteratively 
at the scale of rrig itself). Finally, for the values of the Yukawa couplings of the third 
generation (s)quarks (/ij, hb), we again take the running ones at the scale of the decaying 
particle mass. 

The SUSY breaking mass parameters Mq, M^, Mj, Mi and Mg are taken to be equal 
in all generations. The trilinear breaking parameters of the first and second generation 
are set to zero, i.e. Au^d,e = Ac,s,fi = 0. Furthermore, we simplify and set \Mi\ = M2/2 
(GUT relation), Mq = = M^, M^ = Mg, \At\ = \Ab\ = \Ar\ and if a, = ^a, = ^Ar- 
To make sure that our chosen complex parameter set of the MSSM is not already ruled 
out by the experimental limit on the electric dipole moment of the electron (eEDM), 
we account for the constraint by calculating and thus checking the eEDM automatically 
along the way (see Appendix [F]) . Due to this stringent constraint, we set = and 
focus just on the phase of Aj (the phase of Mi is negligible in our case and thus set to 
zero). 

We study two different sets of processes, considering a) all contributions at full one-loop 
level and b) the two contributions with a gluino in the loop (see Graph 1 and 2 in 
Chapter ED. 

In Fig. 110.11 we show the asymmetry 6'"^ and 6^^ x BR, taking all contributions. We 
vary the input parameter Mq from 500 to 1500 GeV (but show the output parameter 
for better usability) for various tan f3. One can clearly see the threshold of the decay 
ti tg ai ~ 708 GeV, after which the two gluino contributions dominate over 
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Parameter 


Value 


Mao 


800 GeV 


tan f3 


5 


\Mi\ 


M2/2 







M2 


200 GeV 




1000 GeV 


V'm 





Mq 


1000 GeV 


Mn 


Mq 


Mi 


Mq 


Ml 


600 GeV 


Me 


Ml 


1^1 


350 




7r/4 




\M 




VAt 


\Ar\ 


\M 


fAr 


VAt 







A 






Parameter 


Value 


mz 


91.1875 GeV 




80.45 GeV 


cos 6w 


mw/mz 


Gf 


1.16639 X 10-5 GeV-2 




1/127.9 


Us 


9.1045846740771 x 10"^ 


m„ 


53.8 X 10-3 GeV 


rric 


1.5 GeV 




171.4 GeV 


rrid 


53.8 X 10-3 GeV 


rus 


150 X 10-3 GeV 


mi, 


4.2 GeV 


rrie 


0.51099907 X 10-^ GeV 


rriu. 


105.658389 x 10-^ GeV 


rrir 


1777 X 10-3 GeV 



Table 10.1: The input parameters of the MSSM and the SM. 



10.1. A TYPICAL SCENARIO 



45 




Figure 10.1: 6 and 6 x BR as a function of m^-^ for various tan/3, considering all 
contributions (full one-loop). The parameter m^^ is shown for convenience, but the 
parameter actually varied is Mq from 500 to 1500 GeV. 

In Fig. 110.21 we show in detail the dominance of the two gluino contributions over the 
remaining ones by plotting the gluino-to-all ratio. After the threshold at m^^ ~ 708 GeV 
(Mq ~ 719 GeV), the gluino processes account for ~ 98 % of all processes, depending 
on tan/3. The kink at m^^ ~ 1175 GeV (which can be already seen in Fig. 110. ip comes 
from the threshold of ti ^ tx^, where the graph with Xg — t in the stop selfenergy-loop 
and the two graphs with H^{G^) — t — Xs in the vertex correction begin to contribute. 

0.99 

0.98 

0.97 
^ 0.96 
0I0.95 

0.94 

750 1000 1250 1500 

mr [GeV] 




Figure 10.2: Ratio of 6'"^ between the two gluino contributions and all contributions 
(full one-loop) as a function of m^^ for various tan/5. The parameter actually varied is 



Mq from 719 to 1500 GeV. 



Fig. 110.31 shows the comparison of the two gluino contributions by plotting their respec- 
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tive 5'"^ as a function of Mq from 719 to 1500 GeV (m^^ is again shown for convenience). 
Contrary to the expectation that both gluino contributions should dominate due to their 
strong couphng nature, only the contribution with the gluino in the stop selfenergy-loop 
accounts for 6'"^ in a noteworthy manner. The reason why the vertex correction is so 
suppressed lies in the bj — t — coupling of the Graph 2 in Chapter O However, as 
this coupling is embedded in the one-loop vertex correction in a nontrivial way (see 
Eq. (19. 3p ). there exists no simple explanation for this feature. 




700 900 1100 1300 

[GeV] 



1500 



Figure 10.3: Comparison of the two gluino contributions using 

^CP 

function of 

m^^. Label (1) represents the gluino in the selfenergy-loop and (2) the gluino in the 
vertex-correction. The parameter actually varied is Mq from 719 to 1500 GeV. 




Figure 10.4: 6 and 6 x BR as a function of m^^ for various tan/3, considering 
all contributions at full one-loop except the two dominant contributions with a gluino. 
Again, the parameter actually varied is Mq from 500 to 1500 GeV. 

In Fig. llO.4l we plot 6'"^ as well as 6'~^^ x BR as a function of Mq (rrij^ ) , this time taking all 
contributions except the two contributions with a gluino. The kink at m^^ ~ 1175 GeV 
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Figure 10.5: 6 and 6 x BR as a function of m^^ for various (fAt^ considering all 
contributions. As always, the parameter actually varied is Mq from 500 to 1500 GeV. 



is the same feature already seen in Fig. 110.11 and 110. 2[ The kink at m^-^ ~ 750 GeV and 
tan /3 = 35 comes from the closure of the decay channel ti ^ W bi in the graph with 
— bi — W in the vertex correction. Finally, the kink at m^-^ ~ 708 GeV seen in the 
6'"^ X BR plot comes from the branching ratio BR of the decay ti bxt- Since the 
decay channel into is now open, it subtracts a lot from BR. One again, we can see 
that the contributions without a gluino can be neglected, if the decay into a gluino and 
a top is kinematically possible. 

In Fig. 110.51 we present 6'"^ and 6'"^ x BR as a function of Mq {m^^), taking all contri- 
butions for various ipAf Because the complex phase of At (= Ab^r) is the only source of 
CP violation in our chosen scenario, it highly influences 5*"^. Taking ipAt = 0.15 tt, we 
obtain the highest value of 6'-'^ with our parameter set, 5*"^ ~ 24 % at m^-^ = 1500 GeV. 
Fig. 110.61 shows the dependence of S'"^ and S'"^ x BR on tan/? for various Mq, taking 
all contributions. The maximal asymmetry lies around tan 13 ^ 4%. The higher the 
breaking mass parameter Mq (and thus the mass of the decaying particle), the higher 
the asymmetry 6'-'^, because more and more decay channels open up and hence more 
and more processes start to contribute. On the other hand, the more decay channels 
open up, the less is left for the branching ratio BR of the decay ti bxf. 
In Fig. 110.71 we study the dependence of 6'"^ and S*"^ x BR on M2 for various Mq, 
taking all contributions. One can nicely see the closure of the dominating decay channel 
tl —>■ tg, because M2 is related to the gluino mass rrig via GUT relations. The higher 
the mass Mq (wi^J, the later this closure happens. The maximal value of 5*^^ x BR 
with our parameter set is 6^^ x BR^ 3.5 % at M2 ~ 370 GeV and Mq = 1250 GeV. 
Fig. llO.Sl shows the dependence on the absolute value of the trilinear breaking parameter 
I At I for several tan/5. Here we alter our relation of the breaking mass parameters to 
Mu = 1.1 Mq, = 0.8 Mq. Otherwise, we would obtain a non-physical result, because 
of the degeneration of the stop masses (see the main- and off-diagonal elements of the 
sfermion mass matrix in Eq. fl5.ip ) which leads to a singularity in the tj propagator of 
the Graph 1 in Chapter [91 On the other hand, this degeneration would not bother us, 
if we have not used the approximation done in Eq. (17. 8p . If we would calculate of 
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M2 [GeV] M2 [GeV] 



Figure 10.7: 5^^ and 5^^ x BR as a function of M2 for various Mg, taking all contri- 
butions. 
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Eq. (17.41) (and therefore the sum F"*" + T~ in 5'^'^) correctly at full one-loop order using 
renormalization, this singularity becomes harmless. Due to our relation of the breaking 
mass parameters, the dependence of the asymmetry on tan (3 is small. 
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Figure 10.8: S^^ as a function of \At\ for various tan/?, taking all contributions and 
setting Ma = 1.1 Mq, Mj = 0.8 Mq. 

In Fig. 110.91 we plot the dependence of 6'^^ on ip^ty taking (a) all contributions and 
(b) all contributions except the gluino contributions. One can clearly see the periodic 
dependance on ipAf The overall maximum lies at ipAt ~ 0.15 vr with tan/5 = 5. 




(/JAj/tt [rad] ipAt/n[itid] 



Figure 10.9: S^^ as a function of ipAt for various tan/3, taking (a) all contributions and 
(b) all contributions except the gluino contributions. 

Fig. 110.101 shows 6^'^ as a function of |yu| for various tan/3, considering again all contri- 
butions. Because the chargino mass falls below its lower bound, the inner area is ruled 
out by experiment and thus masked grey. For high tan/3 the dependence on be- 
comes rather symmetric, because the low cot /3 in the off-diagonal elements in Eq. (15.41) 
diminishes the dependence on /i. At ~ ±300 GeV we can see the transition of the 
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decaying particle Xi between being wino-like (at higher and higgsino-hke (at lower 
resulting in a different behaviour of the coupling and thus 6'"^. 
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Figure 10.10: S""^ as a function of for various tan/5, considering all contributions. 
The inner area is masked out because the chargino mass falls below its lower bound. 

Fig. 110.111 with S^'^ as a function of (pfj_ demonstrates the stringent constraint on the 
phase of fi, coming from the experimental eEDM-limit. Only a very narrow area is 
allowed, outside this area we have set 6'^^ to zero. As tan/? gets bigger, the allowed 
area reduces even more and therefore we have chosen (y9^ = in our input parameter set 
in the first place. The negligible dependence of 6'"^ on can be explained with the 
very narrow parameter range of y^^. Furthermore, the higher tan/? gets, the lower the 
dependence on /i (due to cot /3 in Eq. (15. 4p ) and thus ip^ becomes. 
For completeness, we also examined the decay of the heavier stop t2 into bxt- Because 
the two stop particles barely mix in our scenario, their masses are very similar (m^^ = 
992.10 GeV and m^^ = 1034.73 GeV). Therefore, the resulting plots are alike, as one 
can see in Fig. 110.121 in comparison with Fig. 110.11 and also in Fig. 110.131 compared with 
Fig. [ml 

Finally, we also investigated the influence on the Yukawa couplings hf and ht taken to 
be running. In our scenario, the difference of the asymmetry 6'"^ taken with running 
Yukawa couplings (tested at two different scales, the mass of the decaying particle and 
1000 GeV) and taken with not running ones is negligible. Only at high values of Mq = 
1500 GeV one obtains a small deviation; not running Yukawa couplings yield a slightly 
higher asymmetry of ~ 5 %. 
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Figure 10.11: (5^^ function of ip^ for various tan f3, considering all contributions. 
6^'^ is set to zero in the regions not allowed by the experimental eEDM-limit. 
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Figure 10.12: 5'-^^ and S'-^^ x BR of the decaying particle ^2 as a function of m^-^ for 
various tan/?, considering all contributions (full one- loop). The parameter m^^ is shown 
for convenience, but the parameter actually varied is Mq from 500 to 1500 GeV. 
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Figure 10.13: 6^^ and 6^^ x BR of the decaying particle t2 as a function of m^^ for 
various tan/5, considering all contributions at full one-loop except the two dominant 
contributions with a gluino. Again, the parameter actually varied is Mq from 500 to 
1500 GeV. 
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10.2 Conclusions 

In this thesis we performed a detailed numerical analysis of the CP violating decay rate 
asymmetry and the quantity S'^^ x BR of the processes ii — > 6 xf and h ^ bxt 
(and their CP transformed counterparts), analyzing the dependence on the parameters 
and phases involved. 

The asymmetry 6^^ rises up to ~ 24%, depending on the point in parameter space. 
The combined quantity 5*"^ x BR reaches up to ~ 3.5 %. 

Finally, we want to comment on the feasibility of measuring this asymmetry at the 
Large Hadron Collider (LHC) at CERN, which will go in operation in summer 2008. 

As the decaying stop is a strong interacting particle, its production cross section will 
be large. Therefore, measurement of our decay rate asymmetry 6'"^ will be possible at 
LHC. However, the precise calculation of the measurability is beyond the scope of this 
thesis, as it involves monte-carlo simulations accounting for all the peculiarities of the 
detector, among other things. 

On the basis of our promising results of the CP violating decay rate asymmetry, we 
suggest that experimenters should search for evidence of these new CP violating asym- 
metries in the MSSM, which can be far beyond the small CP violating effects in the 
SM. These new CP violating sources are not only important in terms of baryogenesis 
but also very interesting for the further understanding of the subatomic world. 
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Appendix A 

Listing of All One-Loop 
Contributions 



Here we specify the complete list of all processes at full one-loop level who can contribute 
to CP violating asymmetries. 




6 6 6 




Xk Xk Xk 




Figure A.l: All vertex contributions who add to the CP violation (Part 1). 
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Figure A. 2: All vertex contributions who add to the CP violation (Part 2). 




Figure A. 3: All stop-selfenergy contributions who add to the CP violation. 
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Figure A. 4: All chargino-selfenergy contributions who add to the CP violation. 



Appendix B 

Passarino— Veltman Integrals 



In this chapter we give the definition of the Passarino- Vehman one-, two-, and three- 
point functions [57] and hst some functions with a special argument set. 



B.l Definitions 

We define the Passarino- Veltman one-, two-, and three-point functions in the convention 
of [SH]. For the general denominators we use the notation 

Then the loop integrals in D = 4 — e dimensions are as follows: 

i/A^. (B.2) 

1 '■.D^ 



and 





Ao{ 


ml) = 


Boipl 


ml, 


mj) = 


b,{pI 


ml, 


mj) = 


Bf,u{pt 


ml, 


mj) = 


Co = 


1 


/a 


c, = 


1 


/A 




1 


/A 



-i/A^. (B.3) 

i / " ^^"^ ^ PUiPiu Bu , (B.5) 



1 



(B.6) 

+ , (B.7) 



•pox)i-p2 

= + Pl^iPluCii + (pi^P2i. + P2fjiPlu)Cl2 + P2^iP2yC22 , (B.8) 

where the C's have (p^, {pi —^2)^,^25 ""^O' "^i' ""^2) their arguments. For further details 
about the coefficient functions, some reductions and relations, some analytical expres- 
sions and some special argument sets see [46] . 
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B.2 Special Argument Set 

Here we list the Passarino-Veltman integrals with a special argument set needed in 
Appendix[Fl All masses of the external particles are set to zero {pf = (pi— ^2)^ = pI = 0) 
and two of the internal particles are the same. 
For the scalar B-function we obtain (see [16] ) 

i?o(0, Ml Ml) = A + log (^^^ (B.9) 

with the real UV-divergence parameter A = - — 7 + log(47r), the Euler-Mascheroni 

Constant 7 ~ 0.577216 and the scale parameter Q. 

Using Feynman parametrization we derive for the scalar C-function 

r (f\ - l-a: + log(x) 

with X = Ml/ Ml and / = (0, 0, 0, M^ Ml, Ml). 
The coefficient functions take the form 

C„(/,.i(A..o.(|l)...i.il±^) ,B.., 

ll-18x + 9x^-2x=^ + 61og(x) 
Cnif) + Cuif) = i2Mo^(i_a;)4 • ^^-13) 

In our special case, the relations Ci ~ C2 and Cn + Cu ~ C22 + Cu hold. 



Appendix C 

Calculation of a Generic Structure 



We show the calculation of a generic structure taking the generic structure I in Sectioning] 
as an example. First, we only take the matrix elements of the three vertices alone: 

^0 = iu{kMPR + g^PLH-h) (C.l) 
= iu{q + v){gfPR + g^PL)v{-q) (C.2) 
M2 = tu{q){g^PR + g^PL)v{-iq + p)). (C.3) 

Then we can write down the matrix element of the selfenergy loop from Aii and Ai2 
using Feynman rules. Because of the closed fermion loop we introduce an additional 
overall factor (—1) (for a detailed derivation see [IS])- We obtain 

where we used the abbreviations 

Do = q^-Ml (C.5) 
D, = [q + pf-Ml (C.6) 
X = TT[{g^Pn + g^PL){^ + M,){g^PR + g^PL){<^ + p + M2)]. (C.7) 

We can modify X to 

X = Tr[^{^ + p){g^g^PR + g^g^PL)+^M2ig^g^PR + g^g^PL) 

+ {i+f)Mi{g'^g^PR + g^g^PL) + M^M^ig^g^PR + g^g^Pi)] 
= 2{q^ + {q.p)){gfg!i + g^g!^) + 2M,M2{gM + gig2) (C.8) 



using the relations P^^ = Pr^l. Pr,lPl,r = and Tt{Pr,l) = 2, Tr{0PR,L) = 2{a.b), 
TtHPr^l) = 0. 

Now we switch from A to D = 4 — e dimensions using the SUSY invariant Dimensional 
Reduction regularization scheme DR resulting in d^g/(27r)^ d^g/(27r)^. Substituting 
this and Eq. fIClSl) into Eq. fIClTl) results in 

M,eu = -I (2^)^ DoD^ + il-P))i9ig2 + gig2) + ^M.M^^g^g!^ + g^g^)] . (C.9) 
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APPENDIX C. CALCULATION OF A GENERIC STRUCTURE 



For the calculation of the loop integrals we use the formalism of Passarino-Veltman 
Integrals introduced in Appendix [Bl Using Eq. flB.3|) . flB.4l) and Eq. flB.5|) reduced with 
the metric tensor g^'^ we get 





1 


i 


(27r)^ 




(47r) 


d^q 


(q.p) 


i 


(27r)^ 


DoDi 


(47r) 


d^q 




i 


(27r)^ 




(4vr) 



11 



(C.IO) 
(C.ll) 
(C.12) 



where the B-functions have (p^,M^,M|) as their arguments. The matrix element then 
becomes 



M 



self 



2 (45oo + ml{B, + Bn)) {g^9^ + 9^9^) 



(47r)2 

+2M,M2(9^92+9i92)Bo 



(47r)2 

+2M,M2(9i92+9i92]Bo 



Ao(M^) + Ao(M|) + (M't + Ml- mi)Bo ) {9^2 + 9i92' 



(C.13) 



where we used the on-shell relation = ml and reduced the term with Bi, Bqq and Bn 
down to Aq and Bq using equations found in [46]. The matrix element of the final graph 
is then derived from Eq. fIC.ll) and flC.lSp 



M 



M 



self 



p2 - Ml 



(4vr) 



-,u{ki){KRPR + KLPL)v{-k2 



(C.14) 



with the form factors 
1 



A 



R 



ml - Mo 



29^ 



A,{Ml) + Ao{Ml) + {Ml + Ml- ml)B, ) {9^92 + 9i92' 



+2M,M2(^9i92+9i92]Bo 
Al = ARiR^L) 

just like in Eq. (lO) . 



(C.15) 
(C.16) 



Appendix D 

Calculation of a Coupling 



As an example for the calculation of couplings we derive the tree level interaction of a 
chargino with a sfermion-fermion pair x~^ff' as specified in Chapter [61 
Since charginos as well as sfermions are mass eigenstates that are formed out of a mixing 
of interaction eigenstates, one has to consider the Lagrangian of interaction eigenstates 
first. The relevant terms are derived from the superpotential and the SUSY gauge 
coupling terms (for a detailed derivation see [IS]). We obtain 

+K{f*^v, + Urrl)^l, + h.c. , (D.2) 
Ceak = t9[iilbL + KrL)X^ + ibltL + fMX- + h.c]. (D.3) 

Then we transform the Lagrangians to mass eigenstates using the relations (15.81 15.131) 
resulting in 

^yL = ht{Rl,i*bL + I^fb,t^j,)v:,xt + h{RU^^^ 

+KiRJ,f:Ur + uAm^Xk + h.c. , (D.4) 

= -9[iRlii*h + KrL)V,\xt 

HRam + R[^>r)Ul^ + h.c] . (D.5) 

Summing up the results we gain 

= {htii,v:,-gi^,,v:,)tibLxt + htR'^v:j),?^^^^ 

+ {hrRl2U;2 - 9RlM,)f:u^Xk + hrU;,i>rTj,Xk + h.c. . (D.6) 

Finally we transform the Weyl spinors to Dirac spinors using the relations flE.9| lE.lOp 
and separate the Lagrangian into a chargino-squark-quark and a chargino-slepton-lepton 
part C^+ff> = C^+qq' + ^^+11'- We get 
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APPENDIX D. CALCULATION OF A COUPLING 



+xtiA^k:pR + Af:PL)tb* + x-,{b^hPr + B^:PL)hi* , (d.?) 

L^^fl, = Vr{j^^,PR + J^tPL)xtn + f{B'i'PR + Bi''PL)Xk^r 

+¥u{^t:PR + KfPLHf* + ru{B',^*PR + BtPL)rK , (D.8) 
with the projection operators Pr^l = (1 ± 7^)/2. 

The abbreviated couphng matrices, which can be further transformed using rriw = 
nit = {htu sin (3) / a/2 and nib = (hbU cos (3) / V2, are 

V2 \niwcos(3 

^ki — "-t-'T^jl ^k2 — : ^ ^k2-^il 5 

V2mv(/ sm p 
Bki = htR%Vk2 - gRa^ki 

= ^(^^V,2R^-V2V,^Ri^ 
V2 \niwsmlJ 

Bi = hbR^l,= J"^' (D.9) 



^/2mw cos P 



A'k^ = KR:2*Uk2-gR:^U. 



kl 



A= ( ^—^Uk2Rj2 - V2UkiRa 

V2 \niw cos p 



K = 0, 

5f = -gVki, 



B'k"- = Ku;,= " m,, (D.io) 

V 2mvy COS p 

just hke in (16.31) and (16.41) . These couphng matrices as weU as the Lagrangians have the 
same convention as in [201. 



Appendix E 

Transformation of Weyl to Dirac 
Spinors 

We define a four component Dirac spinor as follows 

^^fiY (E-l) 



where the two quantities x ^^'^ '0 ai'S two component Weyl spinors, standing for particle 
and anti-particle, respectively. 

We choose the chiral representation so the gamma matrices 7° and 7^ and the projection 
operators Pr^l as well as the adjoint Dirac spinor ^ are: 



7^= , (E.2) 




(E.3) 

^ = vl/^o = (V x) . (E.4) 

With these relations one can now carry out all sorts of transformations including prod- 
ucts of Weyl spinors 

1IJ1X2 = ^iPl*2 = ■ (E.5) 

We define the fermions of the SM, the charginos as well as the neutralinos 

0) ^. = (0^^) r=(^^) r = (rtrt), (E.6) 
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X^=(^^) Xt = ix-,xt) ^'^ = {ll) Xu={xtxl). (E.7) 

^° = fto) M = ixl Xl) . (E.8) 
we obtain for example the following transformations using Eq. (IE. 51) 



^LXk 


= XtPLt 


hxt 


= XkPLh 


tRXt 


= iPLxt 


bnXk 


= hPLXu 


tLXl 


= XlPLt 


hxl 


= XlPLb, 


t\xl 


= iPLxl 




= bPLxl, 



^LXt = XkPLr rlxk = ^PLXk ^rXk = X^Pl^^t , 
rLxl = WLr rlxl = fPLxl = ^^L^r • (E.IO) 



Appendix F 

The Electric Dipole Moment (EDM) 



The additional CP- violating phases in the MSSM are new sources of CP violation beyond 
the SM. From the point of view of baryogenesis, one hopes that these phases are large [6l 
[TJ El [9]. But the experimental limits on electron and neutron electric dipole moments 
(EDMs), 141 < 1-6 X 10-27 [e x cm] [52] and |4| < 2.9 x IQ-^^ [e x cm] [501, place 
constraints on the CP violating phases of the MSSM. Especially the complex phase of /i 
is severely constrained with about (f)^ < 0(10"^) [611[62],[63l|6l| for a typical SUSY mass 
scale of the order of a few hundred GeV. A larger (j)fj_ imposes fine-tuned relationships 
between this phase and other SUSY parameters [Ml EHl l67j . 

In this chapter we calculate the EDM of an electron. We start by giving the relevant 
generic structures and then derive the EDM of a fermion. We apply this result to 
compute the EDM of the electron and point out some numerical issues. Finally we 
use this calculation as an automatic checkup-routine in our calculation of CP violating 
asymmetries so as to make sure our chosen complex parameter set of the MSSM is not 
already ruled out by experiment. 



F.l Generic Structures 



Here we list the two generic structures which are needed for the calculation of the EDM 
of a fermion. Note that we rotated these two generic vertex corrections to fit the generic 
structures found in [IB]. The convention for the momenta and the masses from Chapter [H] 
still apphes but is of course rotated as well. 
The generic matrix element is 

M = j^e^ip)uih){rA, + k'^B.+p'^OPM-h) (F.l) 

with i G {R,L}. The form factors A^, Bi and Cl can be again easily obtained by 

exchanging right- and left-handed couplings {g^ ^ g^) . 

The first generic structure is a scalar-fermion-fermion vertex correction. 
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/2,TO2 




K, Wo 







i{g^PR + g^PL) 



The appropriate form factors are 

^Uf) = - {9o9i9^ (2C00 - Boiml, Ml Mj)) + {9X''h^2'' ' 9^9^i9^Ml)Co 

+ i9^9fh^'' - 9o9^h^'')m,C, - (^7o^^?f /if - g^g^hf'^)m,C2) , (F.2) 
Bkif) = -{9o9^h^''C^+g^g^h^^C2 

+g^g^g^m,{C^, + C12) + 909192^2(022 + C,^)) , (F.3) 

Ckif) = -{9^9^h'^''C^-g^g^h^^C2 

+9o9i92mi{Cu - 0,2) - 9o9i 92^022 - C12)) (F.4) 
with the auxihary functions 

= 9lm, + 9-M, I e {1, 2} ■,j,ke {L, R} (F.5) 
where we do not sum over the index i. The argument set is 

/ = (mo, mi, m2, Mq, Mi, M2, gl^, g^, g^, g^, g^, g^) . (F.6) 
The second generic structure is a fermion-scalar-scalar vertex correction. 



f2,rn2 








2^o(2g + + A;2)^ 



F.2. ELECTRIC DIPOLE MOMENT OF A FERMION 
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The appropriate form factors are 

A'j^if) = 2000909^9^, (F.7) 
B'r'W) = 90 {m.g^g^iCi + Cn + Cu) + ^2(7^(^2 + C22 + 

-M,g^g^iCo + C, + C2)), (F.8) 
C'n'if) = 90 [m.g^g^iCn - Cu) - m2gMiC22 - Cu) 

-MogMiCi - C2)] (F.9) 



with the argument set / = (mo, mi, m2, Mq, Mi, M2, 5-0, 9i, 9i, 92^ 92)- 

F.2 Electric Dipole Moment of a Fermion 

Now we are ready for the calculation of the EDM of a fermion. We start with the matrix 
element in Eq. fIF.ip and apply the Gordon decomposition in order to obtain the a^^- 



terms which describe the spin exchange of the two fermions. Then we use these terms 
to derive the EDM using the relations found in [68] . 
The Gordon decomposition takes the general form 



.0" 



u{p')-i^Piu{p) = —u{p'){p' + pYPMp) + u{p')i—{p' - p)uPMp) (F.IO) 
2m 2m 

with i G {R, L}. Using this relation one can transform the first two term^ in Eq. (IF.ip 
and we obtain 



= 7T^^i^(P)\^i\n u{ki)k^'Piu{k2) +u{ki)i- PuPiu{k2] 

(47r)^ V \2mf 2m f 

^Bi{2mfu{h)^^Piu{k2) - iu{k,)(T^'''p,PMk2))) • (F.ll) 

Taking only the terms with cr^*^ and rearranging a bit yields for the amplitude T [AA"^^" = 
iT) (see [6H]) 

T = tee,{p)^u{k^)a^\a''PR + a'^PL)u{k2) (F.12) 
2mf 

with the coefficients 



= -^^{AR-2mfBR) (F.13) 



e(47r; 



= -^^{AL-2mfBL). (F.14) 



The EDM of a fermion is then calculated by the simple formula^ 



df = -- Im(a-^-a^). (F.15) 

' 4m/ 



^The third term with Ci vanishes in the case of a photon due to the Lorentz condition ef^{p)p^ = 0. 
^On the right-hand-side of Eq. (12) in [68j there is a factor — l/(2my) missing. 
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F.3 Electric Dipole Moment of the Electron (eEDM) 

There are two processes who contribute to the EDM of the electron, one with a sneutrino- 
chargino and one with a selectron-neutrahno in the loop (see Fig. IF.lll . 






Figure F.l: The two processes who contribute to the eEDM. On the left a sneutrino- 
chargino and on the right a selectron-neutralino in the loop. 

The form factors for the sneutrino-chargino contribution have the following arguments 

2 

A]t,L = 5^<L(0,me,me,m,,,m-+,m-+,e,e,i?:^,i?f,i?f*,i?f*) (F.16) 
1=1 

and analogously for Bj^j^. 

The form factors for the selectron-neutralino contribution have the following arguments 
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i=i j=i 



e.,e,D';f,D%D'^,Df) (F.17) 



and analogously for B\i^. 

Because of the huge masses of the particles in the loop compared to the mass of the exter- 
nal electron, one has to carefully avoid numerical instabilities. Especially the Passarino- 
Veltman integrals, where we set the electron mass to zero, need special treatment. As 
now all three external particles are massless, we have to use the special Passarino- 
Veltman integrals derived in Appendix [B] in order to get a numerically stable result. 
After the careful calculation of the form factors we inserted them in Eq. flF.lSP and 
f lRTi]) and calculated the EDM of the electron with Eq. flRlSj) . To compare the result 
with the experimental boundaries given in units of [e x cm] we used the conversion factor 
(me)-^[GeV]-i = 3.861592641934875 x lQ-^^[cm\. 

We checked our numerical result of (using our own coefficients and a^) with the 
result of de obtained by using the coefficients and af^ found in [68]. Finally we incor- 
porated this calculation of the eEDM as an automatic checkup-routine in our calculation 
of CP violating asymmetries in order to verify our complex parameter set of the MSSM. 
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